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1 Introduction 

Almost all of the anti-triplet states A+, S+, S° [ A+(2593), S+(2790), S°(2790)] with 
J^ = I [ I ] ^''^^ containing single heavy charm quark as well as the ^ [ § ] sextet 
Qc, Sc, S^ [Q*, S*, S*] states have been detected in the experiments [1]. Among the S-wave 
bottom baryons, the A^, S^, S^, S^ and Qb states have also been observed. It is expected 
that the LHC not only will open new horizons in the discovery of the excited bottom 
baryons but also it will provide possibility to study properties of heavy baryons as well as 
their electromagnetic, weak and strong decays. 

Such an experimental progress stimulates the theoretical studies on properties of the 
heavy baryons as well as their electromagnetic, weak and strong transitions. The mass 
spectrum of the heavy baryons has been studied using various methods including heavy 
quark effective theory [2], QCD sum rules [3-6] and some other phenomenological models 
[7-12]. Some electromagnetic properties of the heavy baryons and their radiative decays 
have been investigated in different frameworks in [6, 13-24]. The strong decays of the heavy 
baryons have also been in the focus of much attention, theoretically (see for instance [25-27] 
and references therein). 

However, the weak and semileptonic decays of heavy baryons are very important frame- 
works not only in obtaining information about their internal structure, precise calculation 
of the main ingredients of standard model (SM) such as Kabbibbo-Kobayashi-Maskawa 
(CKM) matrix elements and answering to some fundamental questions like nature of the 
CP violation, but also in looking for new physics beyond the SM. The loop level semileptonic 
transitions of the heavy baryons containing single heavy quark to light baryons induced by 
the flavor changing neutral currents (FCNC) are useful tools, for instance, to look for 
the supersymmetric particles, light dark matter, fourth generation of the quarks and ex- 
tra dimensions etc. [28,29]. Some semileptonic decay channels of the heavy baryons have 
been previously investigated in different frameworks (see for instance [30-38] and references 
therein). 

The present work deals with the semileptonic decays of heavy 'E.h(c) and Hw^-, baryons to 
the light S and S baryons. The considered channels are either at loop level described by 
twelve form factors in full QCD or at tree level analyzed by six form factors entering the 
transition matrix elements of the corresponding low energy Hamiltonian. Here, we should 
mention that by the "full QCD " we refer to the QCD theory without any approximation like 
heavy quark effective theory (HQET) so we take the mass of heavy quarks finite. In HQET 
approximation, the number of form factors describing the considered transitions reduce to 
only two form factors [39,40]. The considered processes take place in low energies far from 
the perturbative region, so to calculate the form factors as the main ingredients, we should 
consult some nonperturbative methods. One of the most powerful, applicable and attractive 
nonperturbative methods is QCD sum rules [41,42] and its extension light cone sum rules 
(LCSR) (see for instance [43]). We apply the LCSR method to calculate the corresponding 
form factors in full theory. In this approach, the time ordering multiplication of the most 
general form of the interpolating currents for considered heavy baryons with transition 
currents are expanded in terms of the distribution amplitudes (DA's) of the light S and S 
baryons. Using the obtained form factors, we calculate the decay rate and branching ratio 
for the considered channels. 



The introduction is followed by section 2 which presents the details of the application 
of the LCSR method to find the QCD sum rules for the form factors. Section 3 is devoted 
to the numerical analysis of the form factors as well as evaluation of the decay widths and 
branching fractions. Finally, section 4 encompasses our conclusion. 

2 LCSR for transition form factors 

This section is dedicated to the details of calculations of the form factors. As we previously 
mentioned, the considered transitions can be classified as loop FCNC and tree level decays. 
The loop level transitions include the semileptonic S5 — )■ S/^/^, S;, — )■ E/^/~, Sc — )■ S/^/~ 
S'^ — )■ S/+/^, S^ — )> S/+/^ and H^ — > S/+/^ decays. Considering the quark contents and 
charges of the participant baryons, these channels proceed via FCNC 6 — > s, 6 — > rf or 
c — > M transitions at quark level. The low energy effective Hamiltonian describing the 
above transitions is written as: 

^17/ = ^^^^^^^|C|^^ g7,(l - 75)g/V/ + C,o g7.(l - 75)Q/V75/ 

- 2mQ Cf^^ qta^uQ^il + l5)QlYl\ , (1) 

where Q corresponds to 6 or c quark, Q' represents the t or 6 quark and q denotes the s, 
d or u quark with respect to the transition under consideration. The tree level transitions 
include the channels, Sc — )■ S/z/, Sc — )■ T.lu, E[ — )■ Elu and E'^ — )■ S/z/, which proceed 
via c — )■ s or c — )■ (i depending on the quark contents and charges of the initial and final 
baryons. The effective Hamiltonian representing the considered tree level transitions has 
the following form: 



n'jfj = ^V^e g7^(l - 75)c/r(l - 75)^, (2) 



where q can be either s or d quark, Gp is the Fermi couphng constant, and Vq'q , Vg'q and 
Vqc are elements of the CKM matrix. 

In order to get the amplitudes, we need to sandwich the effective Hamiltonians between 
the initial and final states. Looking at the effective Hamiltonians, we see that we have two 
transition currents, Jjj"'^ = ^7^(1 — 75)Q and Jj['^^ = qi(J fivq'^ {X—lb)Q ■ The matrix elements 
of the transition currents are parameterized in terms of form factors in the following way: 

{B{p) I Jf' I Bq{p + g, s)) = ub{p) [iMQ') + i<y,uq''f2m + q'hm 

- 7m75^i(<5^) - i(^f,ul5q''92{Q^) - q''l593{Q^) ubq{p + q,s) , (3) 



and 



{Bip) I J^'^^ I BQip + g, s)) = ub{p) i,ft m + m^.gV2 m + q'U (Q 



.•'fTfn^\ ,nl^fT( 



+ if^i^gi {Q) + ^(yf.ui^q"gi (Q ) + (i^i^ai (Q ) "bq (p + g, s) , (4) 



where Q^ = — g^, fi, gi, ff and gf are transition form factors, and ubq and ub are spinors 
of the initial and final states. The Bq{p + g, s) stands for particles with momentum p + q 
and spin s. From the explicit expressions of the effective Hamiltonians, it is clear that the 
loop level transitions contain both transition matrix elements having twelve form factors 
while the tree level channels include only the transition current I that corresponds to six 
form factors. 
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Figure 1: Diagrammatic representations of the correlation functions given in Eq. (5). The 
left (right) diagram belongs to the tree (loop) level transitions. The ovals (wavy lines) in 
diagrams stand for the DA's of S or S baryons (external currents). In each graph, the gi 
and 52 are light spectator quarks. 

Our main task in the present work is to calculate the transition form factors. According 
to the philosophy of the QCD sum rules approach, we start with the following correlation 
functions as the main building blocks of the method: 



ni(p,g) 



(fxe"'%B{p) I r{j;'^''(x), 7^^3(0) } I 0) 



Ii';{p,q)=t j d'xe''^%B{p)\T{jy'{ 



a;),J^'?(0)}|0) 



(5) 



where J^'^ is the interpolating current carrying the quantum numbers of the Sq(Sq) 
baryons. The diagrammatic representations of these correlation functions are presented 
in Figure 1. The interpolating currents for the considered baryons have the following gen- 
eral forms (see for instance [44]): 



J-Q = - 



^abc 
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b\ „c 



qr CQ" 75^2 + /3 qr Cl,Q' U- Q"' Cq^ hml + PiQ"' C^^ql ql 



^aT, 



)aT, 



b\ „c 



jHg ^ _,a6c|2('gaT^gy^^gc^2/3(gfC75gyg^+ (gfC7Q''J75g2^ + /3(gfC75g'jg2' 
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Q''^Cqi)^,ql + p{Q'^^C^,qi)ql 



(6) 



where C is the charge conjugation operator, a, h and c are color indices and the light quarks 
gi and g2 are given in Table 1. The /3 is an arbitrary parameter and the value /3 = — 1 
corresponds to the loffe current. 
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Table 1: The light quark contents of the heavy baryons Sq and Sq. 



The correlation functions given above can be calculated in two different ways. From the 
phenomenological or physical side, they are calculated inserting complete sets of hadronic 
states having the same quantum numbers as the chosen interpolating fields. The results of 
this side appear in terms of hadronic degrees of freedom. On the other side, the QCD or 
theoretical side of the correlation functions are calculated in terms of the B baryon DA's 
via operator product expansion (OPE). Then, we match these two different representations 
to relate the hadronic parameters to fundamental QCD degrees of freedom which leads 
to QCD sum rules for the considered form factors. To suppress contribution of the higher 
states and continuum, we apply Borel transformation with respect to the initial momentum 
squared to both sides of the sum rules and use the quark-hadron duality assumption. 

Inserting complete set of hadronic state into correlation functions and isolating the 
contribution of the ground state, we obtain the following representations from physical 
side: 



ni(p,g) 



{B{p) I J^'^ I BQ{p + q,s)){BQ{p + q,s) \ J^<^{0) \ 0) 



Z^ 



m 



Br 



{p + qy 



(7) 



Kip^q) 



y {Bjp) I Jj:''' I BQ{p + q,s)){BQ{p + q,s) \ J^g(O) | 0) ^ 

4^ m| _ (p + g)2 



, (8) 



where the ... stands for the contributions of the higher states and continuum. To proceed, 
besides the transition matrix elements, we need also to know the matrix element {Bq{p + 
q, s) I J^'^ (0) I 0) defined in terms of the residue \bq , 



{BQip + q,s) I J^«(0) I 0) = \BQUBQip + q,s) . 



(9) 



Putting all definitions in Eqs. (7) and (8) and using the completeness relation for Dirac 
particle as 



Yl ubq (p + g, s)ubq {p + g, s) =^+ ^ + ruB^ 



(10) 



we get the following final representations of the correlation functions in physical side: 
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9i m + g's m q, 4l'. 



other structures 



(12) 



where we choose the represented structures to obtain sum rules for the form factors or their 
combinations. Here, we should comment that besides the presented structures, there are 
other structures which one can select to find the form factors. However, our calculations 
show that the selected structures lead to the more reliable results having good convergence 
of sum rules, i.e. in the coefficients of the selected structures, contribution of the higher 
twists is less than those of the lower twists. 

Now, we turn our attention to calculate the QCD sides of the aforesaid correlation 
functions. They are calculated in deep Euclidean region, where —{p + qY ~^ oo. Using the 
explicit expressions of the interpolating currents and contracting out the quark pairs using 
the Wick's theorem, we find 
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Hf/ = — Le"^" / d^xe^''^^ 



2(C)^r,(75)p/3 + (C)«!./3(75)pr, + (C) fir,{l5) p<t> 
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(13) 

2(C75)<^,(/)p/3 



+ (<^75)0/3(/)p^ + (<^75)/3,,(/)p0 



icTi^^q^il - 75) 



s,i-x)p,{o\s;io)s'',ixXMEip)) , 

(14) 



for Eh -)■ SZ+Z" 






2(C)</.r,(75)p/3 + (C)0/3(75)pr, + (C)/3r,(75)p«i 
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for Eb -^ U+l- 
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(C)</./3(75)p,, - (C)/3,7(75)p<^ 



+ /3 



(C^75)<^/3(/) 



P»? 



7^(1 - 75) 5.(-a:)^.(0|5^(0)5^(a:)«^(0)|S(p)) , 



(19) 






(C)0/3(75)pr, - (C)/3r,(75)p<i!, 



+ /3 



(C^75)^/3(/) 



p»7 



- {C-ir))fi^{I)p^ 



i^p.q\^ - 75) 5.(-x)^.(0|s«(0)s^(x)..^(0)|S(p)) , 



(20) 



for S( ^ S/+/-, 



n 



"-^^ / rf^xe"^^^ 
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- (C75)/3r,(/)p0 



(C)</./3(75)p,7 - (C)/3>?(75)p0 



+ /3 



(C^75)<^/3(/) 



P»? 



7^(1 - 75) 56(-:r)^.(0M;;(0)s^(x)rf^/0)|S(p)) , 

- (TO ' 
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n 



// 



'^^^ / rf^xe-^^^ 



(C75)/3r,(/)p0 



(C)^/3(75)pr; - (C)/3^(75)p^ 



+ /3 



(^^75)0/3(1) 



P»? 



io^Xil - 75) 



55(-:r)^.(0M:;(0)4(x)ci^(0)|S(p)) 



(22) 



for s; -^ n+r, 



Ui = -^e"^" I d^xe"^'' 



-I 



(C)</,/3(75)p^ - {C) p^{-i^) p4, 



+ /3 



(^75)<^/3(/) 



pr) 



7^.(1 - 75) \s,{-x)^MK,{^)4{x)dimnp)) 



(23) 






(C)</./3(75)p,7 - (C)/3>?(75)p^ 



+ /3 



(C^75)0/3(/) 
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- (C75)/3r,(/)p0 

for S; ^ S/+/-, 



^a^.g'^ll - 75) 5c(-:r)^.(0|<(0)s^(x)rf^(0)|S(p)) 



(24) 
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for Sr. — )> S/z/ 



C ^ ^OJX, 



n, 



Ve 



e"'"^ / rf^xe"^''^ 
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7,(1 - 75) 5e(-:r);3.(0|<(0)s^(x)rf^(0)|S(p)), 

- cry 



for S^ — )> S/z/, 



n^ = — =e"^^ / rf^xe-*^" 



(C)^/3(75)pr, - (C)/3r,(75)p</, 



+ /3 



(26) 



(C75)0/3(/)p,, - (C75)/3,,(/)p</, 



7,(1 - 75) S,{-x)pMs';{^)4{x)uim'E{p)), 



(27) 



for S^ — )■ S/z/, and 



V2 
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(C75) 0/3 (/),,, - (C75)/3,,(/)p</, 
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(2^ 



for S^ — )> S/z/, where Sq{x) is the heavy quark propagator which is given by [45]: 

l^ + rriQ 



Sq{x) = Sl-{x)-^g,j^^e-^^^ fdv 







{ml - F) 



+ -;^-p^^mG"^'^7. 



m; 



(29) 



and, 



"^Q -4vr2 ^/Z^ ^4^2^2^2(m,V x), {Si)) 

with i^j being the Bessel functions. In Eq. (29), the Sq^ corresponds to the free propaga- 
tion of the heavy quark. The interaction of the heavy quark with the external gluon field 
is represented by the remaining terms. However calculation of these types of interactions 
requires knowledge of the currently unknown four-and five-particle baryonic DA's. The 
contribution of such terms are expected to be small [46-48] , hence, in the present work we 
ignore their contributions. 

To complete the calculations in QCD side, we need also the wave functions of the S 
and S baryons, i.e., e'^^^(0|s^(0)s^(x)M^(0)|S(p)) and e«^^(0|M(rf)^(0)s^(a;)rf^(0)|S(p)). These 
wave functions are expanded in terms of DA's having different twists which are calculated in 
[49] and [50] . For completeness, we present the explicit forms of the wave functions together 
with the DA's in the Appendix. Using the wave functions and heavy quark propagator we 
obtain the correlation functions in QCD side. 

To obtain sum rules for the form factors, we match the coefficients of the same Dirac 
structures from both sides of the correlation functions. We also apply Borel transformation 
and continuum subtraction to suppress the contribution of the higher states and continuum. 
These processes bring us two auxiliary parameters, namely Borel mass parameter M^ and 
continuum threshold sq which we will find the working regions for these quantities in the 
next section. In the meanwhile, we need also the residues Ahq(h' ) whose explicit forms are 
given in [24] . The explicit forms of sum rules for the form factors are very lengthy and we 
do not present their explicit expressions here, but we will give their fit functions in terms 
of g^ in next section. 

3 Numerical Results 

In this section, we numerically analyze the form factors and obtain their behavior in terms 
of q^. Using the fit functions of the form factors, we also calculate the decay rates for all 
considered channels and branching ratios for the channels in which the lifetime of initial 
particle is known. Some input parameters used in the numerical calculations are: m^o = 
(5790.5 ± 2.7) MeV, m^/ = (5790.5 ± 2.7) MeV, m^o = (1314.86 ± 0.20) MeV, m^o = 

(2470.88+IJJ^) MeV, m^!l = (2575.6 ± 3.1) MeV, m^,o = (2577.9 ± 2.9) MeV, m^o = 
(1192.642 ± 0.024) MeV' m^- = (1197.449 ± 0.030) MeV, m^ = (4.7 ± 0.1) GeV, nic = 
(1.27l[J;[!^) GeV, \Vcs\ = 1.023 ± 0.036, iKdl = 0.230 ± 0.011, \VtbVtd^\ = 8.27 x lO'^, 



\VthVtsA = 0.041 , Vbc = (41.2 ± 1.1) X 10-3, Vbu = (3.93 ± 0.36) x lO^^ [1], C^^^ = -0.313, 
Cl^^ = 4.344 and Cio = -4.669 [51]. 

The main input parameters of the LCSR for form factors are the DA's of the S and S 
baryons presented in the Appendix. These DA's contain also four independent parameters. 
These parameters in the case of S baryon are given as [49] : 

/= = (9.9 ± 0.4) X lO--'^ GeV^ Ai = -(2.8 ± 0.1) x 10"^ GeV^ 

A2 = (5.2 ± 0.2) X 10~=^ GeV^ A3 = (1.7 ± 0.1) x 10"^ GeV^ (31) 

and for S baryon, they take the values [50]: 

/s = (9.4 ± 0.4) X lO"-'^ GeV^ Ai = -(2.5 ± 0.1) x 10"^ GeV^ 

A2 = (4.4 ± 0.1) X 10^2 GeV^ A3 = (2.0 ± 0.1) x 10"^ GeV=^. (32) 

The LCSR for form factors contain also three auxiliary parameters. Borel mass param- 
eter M^ and continuum threshold Sq are two of them coming from the Borel transformation 
and continuum subtraction, respectively. The general parameter /3 is the third parameter 
entering the calculations from the general form of the interpolating currents for Bq baryons. 
According to the standard criteria in QCD sum rules, the results of form factors should 
be independent of these auxiliary parameters. Hence, we should look for working regions 
of these parameters such that the dependence of the results on these parameters are weak. 
The working region for the Borel mass parameter is determined requiring that not only the 
higher states and continuum contributions constitute a small percentage of the total dis- 
persion integral but also the series of the light cone expansion with increasing twist should 
converge. This leads to the interval 15 GeV^ < M^ < 30 GeV^ for bottom baryons and 
4 GeV^ < M^ < 10 GeV^ for charmed baryons. The continuum threshold Sq is not totally 
arbitrary but it is related to the energy of the first excited state. Our numerical calculations 
show that in the region (m^Q + 0.3)^ GeV^ < so ^ {''^Bq + 0.7)^ GeV^, the results of the 
form factors exhibit very weak dependency on this parameter. Our numerical calculations 
also lead to the working region —0.6 < cos^ < 0.3 with tan^ = /3 for the general parameter 
/3. As an example, we present the dependence of the form factor /2 for Sf, — )■ 'Ei~^i~ on cos 6 
and M^ in Figures 2 and 3, respectively. From these figures, we see that the form factor /2 
depends weakly on the M^ and Sq compared to the cos^. However, the dependence of the 
/2 on cos 9 in the above mentioned working region is minimal compared to the intervals out 
of the working region. 

Now, we proceed to find the q^ dependence of the form factors in whole physical region, 
i.e. imj < q^ < {tubq — "H^b)^ for loop level and nij < q"^ < {itibq — thbY ^^ ^^^e level 
transitions. However, unfortunately the sum rules for form factors are truncated at some 
points and are not reliable in the whole physical region. This point for instance for the 
Sfe — 7> S/+/^ transition is roughly at q^ = 15 GeV^. To extend the results to whole physical 
region, we look for parametrization of the form factors such that in the reliable region, the 
results obtained from fit parametrization coincide with the sum rules predictions. Using 
the above working regions for the auxiliary parameters as well as other input parameters, 
we find that the form factors are well extrapolated by the fit parametrization, 

M<l'm<l')] = T^^Vt + .i_^.2 - (33) 
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Figure 2: Form factor /2 as a function of cos^ for S^ — )■ Si 
working region of M^ . 
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Figure 3: Form factor /2 as a function of Borel mass parameter M^ for S;, — )■ S/"*"/ decay 
at g^ = 13 GeV"^ and working region of /3. 

The central values for the fit parameters a, b, and mjit as well as values of the form factors 
at g^ = are presented in Tables 7-16. The errors in the values of the form factors at g^ = 
are due to the variation of the auxiliary parameters M^, sq, and /3 in their working regions 
as well as the errors in the other input parameters. To see how the results obtained from 
the fit function coincide well with the sum rules predictions at reliable region, we depict 
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the dependence of the form factors /2 and /^, as examples, on q^ in figures 4 and 5. From 
these figures, we see that the results obtained from the fit parametrization describe well 
the sum rules results in the reliable region. 



c 

n n 


) 


5 1 


D 15 2 




■ ■ ' ■ 






0.2- 






^^^^- 


0.4- 






\ 


0.6- 






- 


0.8- 









10 
q^(GeV^) 



15 



0.0 



-0.2 



-0.4 



-0.6 



-0.8 



20 



Figure 4: Form factor /2 as a function of q^ for S^ — )> Sl^l^ decay at working regions of 
auxiliary parameters. The boxes show the sum rules predictions and the solid line belongs 
to the result obtained from fit parametrization. 



^ 0.8 




q^(GeV^) 



Figure 5: Form factor /J as a function of g^ for S;, — )■ S/^/~ decay at working regions of 
auxiliary parameters. The boxes show the sum rules predictions and the solid line belongs 
to the result obtained from fit parametrization. 

Our last task is to calculate the decay rates and branching ratios of the considered 
channels using the fit functions of the form factors. Considering the amplitudes of the 
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transitions and definitions of tlie transition matrix elements in terms of form factors, the 
differential decay rate for loop level transitions is obtained as [37] : 









81927r5 



0(s) + -A{s) 



(34) 



where s = q^/m\ Gp = 1.17 x 10"^ GeV-^, A = A(l, r,s) with \{a,b,c) = a'^ + b'^ + c^ - 



Arni . 



2ab — 2ac — 2bc and v = \/l 5^ is the lepton velocity. The functions 0(s) and A(s) are 

given as: 



e(s) 



+ 
+ 

+ 



+ 



+ 



32m2m|j^s(l + r - s) {\Dsf + l^gl") 
QAm]m\^{l - r - s) Re[DlE^ + D^El] 
64m|^v^(6m^2 - m|^s)Re[D*^i] 

QAmlm\^y/^(2mBQsRe[DlE^] + (1 - r + s)Re[D*D3 + ElE; 

32m|^(2m? + m|^s)|(l - r + s)mBQV^ Re[A*A2 + 5*^2] 

mB^(l -r-s) Re[A*52 + A^^i] - 2^{Re[AlBi] + m|^sRe[A*52 

l-r)^-s2l}(|Ai|V|5i|^) 

} (1^2!' + 1^21') 

8m|^st;2| - 8mBQSv^Re[D2*^2] + 4(1 - r + s)v^Re[D*D2 + EIE2] 



iTn~Q„ <; 4m Al + r — s 



im^ i Am 



'^rn^Bn > 4m^(l + r — s 



A + (1 + r- s)s 



"^Iq^ 



'l-r?-s' 



'l-r?-s^ 



- 4(1 - r - s) Re[D*E2 + D^Ei] + m^. 



\2 2 

r) — s 



{\D,\' + \E,f)] , (35) 



A (s) = -Sm^'syX (lAil' + IfiiT + l^il + l^il) 

' \E2\ 



+ 8m%^sv^x(^\A2f + \B2f + |/^2 



where r 



m.\/m.\>^ and 



'Bq 



1 



^1 = - (/r + al) {-2mQCr) + (/i - gi] 

A2 = A (1 ^ 2) , 

A3 = Ai (1 ^ 3) , 

El = Ai (^1 -^ -gi; gl -^ -gj) , 

B2 = 5i (1 ^ 2) , 

E3 = 5i (1 ^ 3) , 

/^i = (/i-^i)Cio, 



'^9 



(36) 
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D2 = Di (1 ^ 2) , 

D3 = /^i (1 -^ 3) , 

El = Di {gi -^ -gi) , 

E2 = El (1 ^ 2) , 

E, = El (1 ^ 3) . (37) 

Integrating the differential decay rate over s in whole physical region, 4m'j/m'^^ < s < 
(1 — v^)^5 0116 can obtain the total decay rate. 

For the tree level transitions, the formula for the decay width is given by [52,53]: 

-- ^^ 'V;s(d)P / dq' (1 - m^/qy ^{a^ - q^W ' Q') N{q') (38) 



-- mf 

where 

Niq') = Fl{q'){5\Aq^-m]) + 2a^5\l + 2myq^)-{a^ + 2q'){2q^ + m])) 

+ F'^iq'W - q')i2a' + q'){2q' + mf)^,, + 3F|(g^)mf (a^ - g^)gV^^ 



(') 



+ 6Fi(g2)F2(g2)(52 - q^)(2q' + mj)a/m^,,, - m{q')F,{q^Wi{a^ - q^)6/m^,,, 

+ G'2(g2)(a2(4^2 _ ^2) ^ 2^2^2(1 + 2m2/g2) _ (52 ^ 2g2)(2g2 + m'l)) 

+ G^(g2)(a2 - g^)(2<52 + q'){2q' + mf)/m^„ + 3G^(g^)mf (5^ - g^)gV^l,') 

- 6Gi(g2)G2(g^)(a2 - ^2^(2^2 ^ ^2^^/^^^,^ ^ 6Gi(g2)G3(g')mf (5^ - g2)a/m^n , 

(39) 

with Fi(g2) = /,(g2)^ ^^(^2) ^ ^^,,/2(g2)^ ^3(^2) ^ ^^„/3(g2)^ ^^(^2) ^ ^^(^2)^ ^^(^2) ^ 

"^=.(')fl'2(Q'^), ^3(5^) = 'rn„(i)g^{q^), a = m„(i) + m^, 5 = m„(') — m^ and mi is the lepton's 

^c ^c ^c ^c 

mass. The numerical results of decay width for considered channels are presented in Ta- 
ble 17. Finally, for the channels which we know the lifetime of the initial particles [1], we 
calculate the branching ratios as presented in Table 18. The orders of branching fractions 
for most of the channels presented in Table 18 show that these channels are accessible at 
LHC. 



4 Conclusion 

In the present study, we have considered various loop level and tree level semileptonic 
decays of heavy '^'u^\ and S^c) baryons to the light S and S baryons in the framework of 
the light cone QCD sum rules. The most general form of the interpolating currents for 
the considered heavy baryons as well as the recently available distribution amplitudes of 
the S and S baryons have been used to calculate twelve form factors for loop level and six 
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form factors for tree level transitions in full theory of QCD. Using the sum rules for the 
form factors, then, we have evaluated the decay rates of the related transitions. For those 
transitions with known lifetime, we have also calculated their branching fractions. The 
orders of branching fractions for tree level S^ — )> S l^ui and Sc — )> S l^ui (with I = e or fi) 
as well as rare loop level Sf, — )> S l^l^ and S^ — ?■ S /+/^ (with / = e or /i or r) transitions 
show that these channels can be detected at LHC The similar baryonic A5 — )■ A fi~^fi~ has 
been observed very recently by CDF Collaboration [54] and they reported the branching 
ratio of [1.73±0.42(stat) ±0.55(syst)] x 10~® which is in good consistency with our previous 
work [37]. Any measurement on the considered channels in the present work and comparison 
of the obtained data with our results can help us understand better the internal structures 
of the considered heavy baryons as well as obtain useful information about the distribution 
amplitudes of the S and S baryons. Such comparison in FCNC channels can help us also 
in the course of searching for new physics effects beyond the SM. 

5 Acknowledgment 

Two of the authors (K. A. and H. S.) would like to thank TUBITAK, for their partial 
financial support provided under the project No.llOT284. 



14 



References 

[1] K. Nakamura et al, J. Phys. G 37 , (2010) 075021. 

[2] A. G. Grozin, O. I. Yakovlev, Phys. Lett. B 285 (1992) 254; Y. B. Dai, C. S. Huang, 
C. Liu, C. D. Lu, Phys. Lett. B 371 (1996) 99; D. W. Wang, M. Q. Huang, C. Z. Li, 
Phys. Rev. D 65 (2002) 094036. 

[3] D. W. Wang, M. Q. Huang, Phys. Rev. D 68 (2003) 034019; Z. G. Wang, Eur. Phys. 
J. C 54 (2008) 231; F. O. Duraes, M. Nielsen, Phys. Lett. B 658 (2007) 40; X. Liu, H. 
X. Chen, Y. R. Liu, A. Hosaka, S. L. Zhu, Phys. Rev. D 77 (2008) 014031. 

[4] E. V. Shuryak, Nucl. Phys. B 198 (1982) 83. 

[5] Zhi-Gang Wang, Phys. Lett. B 685 (2010) 59. 

[6] T. M. Ahev, K. Azizi, A. Ozpineci, Nucl. Phys. B 808, (2009) 137. 

[7] M. Karhner, H. J. Lipkin, Phys. Lett. B 660 (2008) 539. 

[8] M. Karliner, B. Kereu-Zura, H. J. Lipkin, J. L. Rosner, arXiv:0706.2163 [hep-ph]. 

[9] D. Ebert, R. N. Faustov, V. O. Galkin, Phys. Rev. D 72 (2005) 034026. 
[10] N. Mathur, R. Lewis , R. M. Woloshyn, Phys. Rev. D 66 (2002) 014502. 
[11] J. L. Rosner, Phys. Rev. D 75 (2007) 013009. 

[12] Qiao-Yan Zhao, Dan Zhang, Qiu-Yang Zhang, Chin. Phys. Lett. 28, (2011) 071201. 
[13] A. L. Choudhury, V. Joshi, Phys. Rev. D 13 (1976) 3115. 
[14] D. B. Lichtenberg, Phys. Rev. D 15 (1977) 345. 
[15] R. J. Johnson, M. Shah-Jahan, Phys. Rev. D 15 (1977) 1400. 
[16] S. K. Bose, L. P. Singh, Phys. Rev. D 22 (1980) 773. 
[17] L. Y. Glozman, D. O. Riska, Nucl. Phys. A 603 (1996) 326. 
[18] B. Julia-Diaz, D. O. Riska, Nucl. Phys. A 739 (2004) 69. 
[19] S. Scholl, H. Weigel, Nucl. Phys. A 735 (2004) 163. 
[20] Y. S. Oh, B. Y. Park, Mod. Phys. Lett. A 11 (1996) 653. 
[21] A. Faessler et. al, Phys. Rev. D 73 (2006) 094013. 
[22] S. L. Zhu, W. Y. Hwang, Z. S. Yang, Phys. Rev. D 56 (1997) 7273. 
[23] T. M. Aliev, K. Azizi, A. Ozpineci, Phys. Rev. D 77 (2008) 114006. 
[24] T. M. Aliev, K. Azizi, A. Ozpineci, Phys. Rev. D 79, (2009) 056005. 

15 



[25] T. M. Aliev, K. Azizi, M. Savci, Phys. Lett. B 696 (2011) 220. 

[26] T. M. Aliev, K. Azizi, M. Savci, Nucl. Phys. A 852 (2011) 141 

[27] T. M. Aliev, K. Azizi, M. Savci, V. S. Zamiralov, Phys. Rev. D 83, (2011) 096007. 

[28] G. Buchalla, G. Hiller and G. Isidori, Phys. Rev. D 63 (2000) 014015. 

[29] C. Bird, P. Jackson, R. Kowalewski, M. Pospelov, Phys. Rev. Lett. 93 (2004) 201803. 

[30] D. Ebert, R. N. Faustov, and V. O. Galkin, Phys. Rev. D 73, (2006) 094002. 

[31] Ming-Qiu Huang, Dao-Wei Wang, Phys. Rev. D 69, (2004) 094003; arXiv:0608170v2 
(hep-ph). 

[32] C. Albertus, E. Hernandez, and J. Nieves, Phys. Rev. D 71, (2005) 014012. 

[33] Ruben Flores-Mendieta, J.J. Torres, M. Neri, A. Martinez, A. Garcia, Phys. Rev. D 
71, (2005) 034023. 

[34] Muslema Pervin, Winston Roberts, and Simon Capstick, Phys. Rev. C 72, (2005) 
035201. 

[35] K. Azizi, M. Bayar, A. Ozpineci, Y. Sarac, Phys. Rev. D 80, (2009) 036007. 

[36] K. Azizi, M. Bayar, Y. Sarac, H. Sundu, Phys. Rev. D 80, (2009) 096007. 

[37] T. M. Aliev, K. Azizi, M. Savci, Phys. Rev. D 81, (2010) 056006. 

[38] K. Azizi, M. Bayar, M. T. Zeyrek, J. Phys. G: Nucl. Part. Phys. 37, (2010) 085002. 

[39] T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. B 355 (1991) 38. 

[40] T. M. Aliev, A. Ozpineci, M. Savci, Phys. Rev. D 65 (2002) 115002. 

[41] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[42] M. A. Shifman, A. I. Vainstein, V. L Zakharov, Nucl. Phys. B 147, 448 (1979). 

[43] I. I. Balitsky, V. M. Braun, A. V. Kolesnichenko, Nucl. Phys. B 312 (1989) 509. 

[44] E. Bagan, M. Chabab, H. G. Dosch, S. Narison, Phys. Lett. B 278, (1992) 369. 

[45] I. I. Balitsky, V. M. Braun, Nucl. Phys. B 311 (1989) 541. 

[46] V. M. Braun, A. Lenz, M. Wittmann, Phys. Rev. D 73 (2006) 094019. 

[47] V. M. Braun, A. Lenz, N. Mahnke, E. Stein, Phys. Rev. D 65 (2002) 074011; A. Lenz, 
M. Wittmann and E. Stein, Phys. Lett. B 581, 199 (2004). 

[48] V. Braun, R. J. Fries, N. Mahnke and E. Stein, Nucl. Phys. B 589 (2000) 381. 

[49] Yong-Lu Liu, Ming-Qiu Huang, Phys. Rev. D 80, (2009) 055015. 

16 



[50] Y. L. Liu, M. Q. Huang, Nucl. Phys. A 821, (2009) 80. 

[51] A. J. Buras, M. Muenz Phys. Rev. D52 (1995) 186. 

[52] A. Faessler, T. Gutsche, B. R. Holstein, M. A. Ivanov, J. G. Korner, V. E. Lyubovitskij, 
Phys. Rev. D 78, (2008) 094005. 

[53] H. Pietschmann, Acta Phys. Austriaca Suppl. 12, (1974) 1. 

[54] T. Aaltonen et al (CDF Collaboration), Phys. Rev. Lett. 107, 201802 (2011). 



17 



Appendix A 

In this Appendix, we present the general decomposition of the wave functions of the baryons 
in final states, i.e., e"''^(0|gi^(0)g2e(a;)g3^(0)|5(p)) and DA's of the E and S baryons [49,50]: 



me'"''qiUaix)q2Ua2x)q3':,iasx)\Bip)) 
+ VimBii5C)a^B^ + V2ml{^5C)a(si^By, + (Vl + 



2^2 



X m 



By-iM 



K){^cu{i^B) 



+ V2mBi^C)^^{^-f5B)^ + VsmBi-f^,C)^^i-f''-f5B)j + Vml{^C)^p{-i^B). 



+ 



2 2 

xmt: 



AT){HC)a^B^ + A2mB{HC)ap{^B)^ + ^3"^B(7M75C)a/3(7^5)' 



A4m%{:jt-f5C)afiB^ + Ar>m\{-i ^-i^C)ap{i(y^'' XyB)^ + AQm%{:fe-ir^C)ap{i:B). 



2™2 



+ (ri + 



xm 



Bn-M 



TrWia^,C)^p{-1>'-1^B)^ + r2mB(x>'^2a^,C)„;3(75fi)^ 



+ TzmB{(y^uC)af}{cr^''l^B)^ + %mB{p''(r^^C)af3{<y^^Xp-f5B)^ 
+ r7m|(a^,C)«;3(a^'^;f755)^ + r8m|(xV^,C)«^(a^^x,755)- 



(Al) 



The calligraphic functions in the above expression have not definite twists but they can 
be written in terms of the B distribution amplitudes (DA's) with definite and increasing 
twists via the scalar product px. The relationship between the calligraphic functions ap- 
pearing in the above equation and scalar, pseudo-scalar, vector, axial vector and tensor 
DA's for B baryon are given in Tables 2, 3, 4, 5 and 6, respectively. 



Si = Si 



2pxS2 = Si - S2 



Table 2: Relations between the calligraphic functions and B scalar DA's. 



Vi = Pi 



2pxV2 = Pl-P2 



Table 3: Relations between the calligraphic functions and B pseudo-scalar DA's. 

Every distribution amplitude, F= 81^2, Pi,2, Vi^e, ^1^6) ^i^s can be represented as: 

F{aipx) = / dxidx2dx^5{xi + X2 + x-i- l)e~*P''(^?=i ''i''^^F{xi) ■ {A.2) 

where, Xi with i = 1, 2 or 3 are longitudinal momentum fractions carried by the partici- 
pating quarks. 



2pa;V2 = Vi-V2-V^ 



2V, = V. 



ApxVi 



-2\/i + V3 + V4 + 2V5 



4pxV5 = V; - \/3 



AipxfVe = -Vi + V2 + Vs + V4 + V,-Vg 



Table 4: Relations between the calligraphic functions and B vector DA's. 





Ai 


= Ai 




2pxA2 


= — , 


A1 + A2- 


As 




2A3 


= As 




ApxAi = - 


-2Ai 


-A3-A4 


+ 2A, 


ApxA5 = 


--A3-A4 




Aipx^Ae = Ai 


-A2 


+ A3 + A. 


i-A^ + Ae 



Table 5: Relations between the calligraphic functions and B axial vector DA's. 



The explicit expressions for the DA's of the B baryon up to twists six are given as 
[49,50]: 
Twist-3 distribution amplitudes: 

Vi{Xi) = 120X1X2X303, Ai(Xi) = 0, 

Ti(Xi) = 120X1X2X303°. (40) 

Twist-4 distribution amplitudes: 



S'i(Xj) = 6(x2-Xi)x3(^° + ^4°), 



24xiX204 , 
12X3(1 -X3)^4°, 



24XiX204° , 



V2{X,) 

Vsix,) 

T2{Xi) 

T-j{xi) = 6x3(1 -X3)(e;°-e4°)- 

Twist-5 distribution amplitudes: 

^2(0:.) = \{x,-x2m+^,'), 

V^ix,) = 3(1-X3)V'5, 
V5(Xi) = 6X305, 

3 



--(:ri+X2)(e5'^ + ei 



Tiixi) - 

Tsixi) ■ 
Twist-6 distribution amplitudes: 



'' ' 'I) , 



^i^i + X2M' - e,) ■ 



Ve{x^) 
Tfifxi) 



k'O 



Pl(x,)=6(x2-Xi)x3(e4°-el°), 

A2{xi) = 0, 

^3(3;*) = -12X3(xi - X2)^4 , 
T3(X,)= 6X3(1 -X3)(e4° + e4°), 



P2ix^) = ^{X, - X2M - ^',') , 

A4(Xi) = 3(xi -X2)V'5 , 

A5{x,) = 0, 

T5{Xi) = 6X305°, 



Aeixi) = 0, 



(41) 



(42) 



(43) 
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ri = Ti 




2pxT2 


= Ti + T2 


-2T3 


2% = Tr 




2pxT4 


= Ti - T2 


-2T7 




2px% -- 


= -Ti + T5 + 2Ts 


Kpxf% = 


--2T2- 


2T3 - 2T4 + 2T5 + 2T7 + 2T8 




ApxTr = Tr- 


Tg 


^pxfTs 


= -Ti + T2 + Tg - 


- Tg + 2T7 + 2Ts 



Table 6: Relations between the calligraphic functions and B tensor DA's. 



where, 



4 - 


= 0^ = /b, 


r, = r, = l{fB-\i 


02 = 


= 4 = lifB + \l), 


0? = 0;° = -e, = i( 


01° = 


= e4 = ^(8A3-3A2), 
6 


0'5° = -a' = 1^2 , 


ef = 


= i(i2A3-5A2). 





(4A3 - A2) 



(44) 
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a 


b 


rrifit 


q^ = 


/l 


0.166 


-0.024 


5.35 


0.142 ±0.036 


/2 


0.028 


-0.048 


5.31 


-0.020 ±0.005 


/3 


-0.004 


-0.006 


5.37 


-0.010 ±0.002 


9i 


0.106 


0.054 


5.24 


0.160 ±0.042 


92 


-0.005 


-0.004 


5.28 


-0.009 ±0.002 


93 


0.003 


-0.006 


4.70 


-0.003 ±0.001 


fl 


0.127 


-0.129 


5.10 


-0.0020 ± 0.0005 


fl 


0.072 


0.085 


5.40 


0.157 ±0.041 


fl 


-0.003 


0.049 


5.23 


0.046 ±0.011 


9l 


0.288 


-0.312 


4.80 


-0.024 ±0.006 


9l 


0.036 


0.119 


4.70 


0.155 ±0.040 


9l 


0.024 


-0.095 


5.33 


-0.071 ±0.018 



Table 7: Parameters appearing in the fit function of the form factors and the values of the 
form factors at q^ = for S^ — )■ Si'^i" . 
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a 


b 


rrifit 


q^ = 


/l 


0.035 


0.011 


5.16 


0.046 ±0.011 


/2 


0.027 


-0.060 


5.32 


-0.033 ±0.008 


/3 


0.086 


-0.110 


5.38 


-0.024 ±0.006 


9i 


0.047 


0.020 


5.34 


0.067 ±0.017 


92 


-0.003 


-0.021 


5.25 


-0.024 ±0.006 


93 


-0.003 


-0.024 


5.39 


-0.027 ±0.006 


fl 


0.045 


-0.047 


5.29 


-0.0020 ± 0.0005 


fl 


0.034 


0.015 


5.25 


0.049 ±0.012 


fl 


-0.145 


0.168 


5.17 


0.023 ±0.006 


9l 


0.006 


-0.012 


4.67 


-0.006 ±0.001 


9l 


-0.041 


0.054 


5.38 


0.013 ±0.003 


9l 


0.049 


-0.071 


5.36 


-0.022 ±0.005 



Table 8: Parameters appearing in the fit function of the form factors and the values of the 
form factors at q^ = for S^ — )■ S£"^£~. 
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a 


b 


rrifit 


g^ = 


A 


0.526 


-0.116 


1.53 


0.409 ±0.106 


/2 


-0.550 


0.026 


1.58 


-0.524 ±0.136 


/3 


-0.204 


-0.582 


1.57 


-0.786 ± 0.204 


9i 


0.183 


0.154 


1.55 


0.337 ±0.088 


92 


-0.431 


0.045 


1.63 


-0.386 ±0.100 


93 


-0.190 


-0.285 


1.63 


-0.475 ±0.123 


fl 


0.042 


-0.048 


1.56 


-0.006 ±0.001 


fl 


0.585 


-0.125 


1.52 


0.460 ±0.120 


fl 


-0.449 


1.127 


1.59 


0.678 ±0.176 


9l 


0.058 


-0.062 


1.58 


-0.004 ±0.001 


9l 


0.730 


-0.201 


1.57 


0.529 ±0.260 


9l 


-0.531 


-0.148 


1.61 


-0.679 ±0.176 



Table 9: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for Sc — )■ S£+£~. 
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a 


b 


mjit 


g2 = 


/l 


0.092 


-0.003 


5.30 


0.089 ±0.022 


/2 


-0.010 


-0.021 


5.32 


-0.031 ±0.007 


/3 


0.015 


-0.058 


5.73 


-0.043 ±0.010 


91 


-0.421 


0.477 


5.20 


0.056 ±0.014 


92 


-0.012 


-0.008 


5.10 


-0.020 ±0.005 


93 


-0.035 


0.001 


5.00 


-0.034 ±0.008 


fl 


-1.126 


1.124 


5.40 


-0.0020 ± 0.0005 


fl 


0.028 


0.081 


4.80 


0.109 ±0.028 


fl 


0.035 


0.132 


5.26 


0.167 ±0.043 


9l 


0.645 


-0.645 


5.40 


0.000 ±0.000 


T 
92 


0.022 


0.002 


4.80 


0.024 ± 0.006 


T 

9i 


-0.210 


-0.058 


5.32 


-0.268 ± 0.070 



Table 10: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for 'E', — )■ 'E.(i^(i~ . 
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a 


b 


mfit 


g2 = 


A 


0.034 


0.056 


5.13 


0.090 ±0.022 


/2 


0.046 


-0.104 


5.34 


-0.058 ±0.014 


/3 


0.055 


-0.104 


5.27 


-0.049 ±0.012 


91 


-0.237 


0.275 


5.36 


0.038 ±0.010 


92 


0.008 


-0.049 


5.34 


-0.041 ±0.011 


93 


-0.006 


-0.039 


5.31 


-0.045 ±0.011 


fl 


0.458 


-0.458 


5.15 


0.000 ±0.000 


fl 


-0.541 


0.679 


5.35 


0.138 ±0.036 


fl 


-0.281 


0.494 


5.38 


0.213 ±0.055 


9l 


0.722 


-0.725 


5.08 


-0.003 ±0.001 


T 
92 


-0.106 


0.191 


5.28 


0.085 ±0.021 


T 

9i 


0.025 


-0.327 


5.32 


-0.302 ± 0.078 



Table 11: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for S. — )■ S£+£~. 
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a 


b 


rrifit 


g2 = 


/l 


-0.564 


0.640 


1.52 


0.076 ±0.019 


/2 


-0.426 


-0.258 


1.55 


-0.684 ±0.178 


fs 


-0.642 


-0.297 


1.58 


-0.939 ± 0.244 


91 


-0.092 


0.212 


1.62 


0.120 ±0.031 


92 


-0.265 


-0.081 


1.60 


-0.346 ± 0.090 


93 


0.238 


-0.349 


1.55 


-0.111 ±0.029 


fl 


0.272 


-0.293 


1.60 


-0.021 ±0.005 


/J 


0.432 


0.112 


1.53 


0.544 ±0.141 


fl 


-0.433 


0.605 


1.62 


0.172 ±0.045 


9l 


0.258 


-0.265 


1.72 


-0.007 ±0.002 


9l 


0.401 


-0.013 


1.50 


0.388 ±0.101 


9l 


0.153 


-0.510 


1.63 


-0.357 ±0.093 



Table 12: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for S^ — )> S£+£~. 





a 


b 


rrifit 


g2 = 


/l 


-0.4142 


0.608 


1.52 


0.194 ±0.050 


/2 


-0.320 


-0.036 


1.60 


-0.356 ±0.092 


h 


1.068 


-1.530 


1.55 


-0.462 ±0.120 


91 


-0.624 


0.935 


1.58 


0.311 ±0.081 


92 


0.010 


-0.161 


1.63 


-0.151 ±0.038 


93 


1.398 


-1.710 


1.61 


-0.312 ±0.081 



Table 13: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for Sc — )■ S£z/. 
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a 


b 


rrifit 


g2 = 


A 


0.528 


-0.119 


1.52 


0.409 ±0.106 


/2 


-0.564 


0.006 


1.58 


-0.558 ±0.145 


/3 


-0.598 


-0.133 


1.57 


-0.731 ±0.190 


91 


0.448 


-0.031 


1.55 


0.417 ±0.104 


92 


-0.493 


0.096 


1.61 


-0.397 ±0.099 


93 


-0.295 


-0.329 


1.63 


-0.624 ±0.150 



Table 14: Parameters appearing in the fit function of the form factors and the values of the 
form factors at q^ = for Sg -^ S£z/. 





a 


b 


rrifit 


g2 = 


A 


-1.498 


2.075 


1.60 


0.577 ±0.150 


/2 


-0.359 


-0.142 


1.66 


-0.501 ±0.130 


/3 


-0.760 


0.082 


1.70 


-0.678 ±0.176 


91 


0.159 


0.292 


1.57 


0.451 ±0.113 


92 


-0.317 


-0.024 


1.62 


-0.341 ±0.089 


93 


0.976 


-1.218 


1.62 


-0.242 ±0.061 



Table 15: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for S^ — > Eiu. 





a 


b 


rrifit 


g2 = 


/l 


-0.564 


0.640 


1.52 


0.076 ±0.020 


/2 


-0.226 


-0.427 


1.55 


-0.653 ±0.169 


/3 


-1.007 


0.112 


1.58 


-0.895 ±0.232 


91 


-0.017 


0.054 


1.62 


0.037 ±0.009 


92 


-0.265 


-0.081 


1.60 


-0.346 ± 0.089 


93 


0.238 


-0.349 


1.55 


-0.111 ±0.028 



Table 16: Parameters appearing in the fit function of the form factors and the values of the 
form factors at g^ = for S^ — )■ S£z/. 
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1(^6^) 


Eb -^ Ee+e' 


(9.941 ± 2.982) x IQ-^^ 


Eb -^ Efi+ij,- 


(9.872 ± 3.455) x lO-^^ 


Eb -^ Et+t- 


(1.611 ±0.483) X 10-19 


Eb — )■ Se+e^ 


(6.359 ± 2.225) x lO'^" 


Eb -)■ S;U+;U- 


(6.194 ±2.167) X 10-20 


Eb -^ T.T+T- 


(4.338 ±1.518) X 10-20 


Ec -^ Se+e" 


(2.162 ±0.757) X lO-^^ 


He -^ S/i+/i" 


(2.152 ±0.753) X lO-^^ 


He -^ Se+z/e 


(4.264 ±1.49) X 10-13 


He -^ H/i+Z/^ 


(4.202 ±1.26) X 10-13 


He -)> Se+Z/e 


(2.204 ±0.771) X 10-1^ 


He -> T^fl+U^ 


(2.183 ±0.764) X 10-^^ 


% ^ He+e- 


(2.411 ±0.843) X 10-1^ 


H', -> H/i+/i- 


(2.407 ±0.842) x 10-^^ 


H; -> Hr+r- 


(1.199 ±0.419) X 10-1^ 


K -^ s//+z/^ 


(1.009 ±0.303) X 10-12 


H', ^ He+z.e 


(1.109 ±0.388) X 10-12 


H; ^ Se+z/e 


(8.425 ± 2.948) x lO-i^ 


H:^ ^ S/i+z/^ 


(8.340 ±2.919) X lO-i^ 


H;^ ^ Se+e- 


(1.718 ±0.515) X 10-24 


H', ^ S/i+/i- 


(1.666 ± 0.499) X 10-2^ 


H; ^ Ee+e- 


(3.815 ±1.335) X 10-19 


H; -> S/i+/i- 


(3.813 ±1.334) X 10-19 


H', ^ Sr+r- 


(2.783 ± 0.974) x lO-i^ 



Table 17: The values of the decay rates in full theory for different leptons. 
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BR 


Eb -^ Ee+e- 


(2.25 ±0.78) X 10-6 


Eb -> Efi+/j' 


(2.23 ±0.67) X 10-6 


Eb -^ Et+t- 


(0.36 ±0.11) X 10-6 


Eb -^ Se+e- 


(1.44 ±0.50) X 10-^ 


Eb -^ 1:12+ jj,- 


(1.40 ± 0.49) X 10-'' 


Eb -^ T.T+T- 


(0.98 ±0.29) X 10-^ 


Se -^ Se+e- 


(3.68 ±1.29) X 10-1^ 


Ec -> T.ji^ii' 


(3.66 ±1.28) X 10-^^ 


Sc -^ Se+z/e 


(7.26 ±2.54) X 10-2 


Sc -^ S/i+z/^ 


(7.15 ±2.50) X 10-2 


Ef. — )■ Ee+fe 


(1.48 ±0.52) X 10-2 


He ^ S;U+z/^ 


(1.47 ±0.51) X 10-2 



Table 18: The values of the branching ratios in full theory for different leptons. 
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